Multiscale wireless communication

ABSTRACT

Multiscale communication, wherein time-frequency conditioning of the transmitted signal is used, is examined with respect to typical wireless communications systems, in which the modulating wavelet is matched to a specific channel condition. Compactly supported wavelet bases are employed, and based on their parameterization a wavelet is selected that best characterizes the wireless channel conditions.

FIELD OF THE INVENTION

The field of the invention is that of wireless communications, in particular that of wavepacket systems that employ separation in both the time and frequency domains.

BACKGROUND OF THE INVENTION

Wireless communications involves several forms of signal modulation prior to transmission through a mobile channel. Some examples of the type of processing involved in wireless modulation include temporal processing (e.g. spread spectrum), spectral processing (e.g. orthogonal frequency division multiplexing, i.e. OFDM) and spatial processing (e.g. space-time coding). This kind of processing occurs over a single dimension (e.g. time, frequency, or spatial), and is sometimes referred to as single-scale modulation. Usually the type of processing involved is oftentimes selected based on the type of mobile channel conditions experienced. For instance, a mobile channel that provides no diversity (e.g. single-path fading) is sometimes addressed using spatial processing. On the other hand, a mobile channel where multipath is experienced (e.g. frequency-selective fading) may be better handled using spectral modulation such as OFDM.

Multiscale modulation involves processing the signal over two dimensions, namely time and frequency. Therefore, the output of a multiscale modulator is indexed by both a temporal range and frequency bin. This type of signal conditioning has the potential to match to two dimensions of the wireless channel rather than one.

Multiscale modulation can be visualized using a time-frequency tiling diagram. A sample tiling, derived from (1), is depicted in FIG. 1. The time-frequency tiling of a waveform composed of a sinusoid at frequency f₀ and impulse t₀ in this example results in energy in all subbands in the time-frequency domain intersecting both f₀ and t₀. A time-frequency representation of a signal can be obtained using a wavelet decomposition.

Previous work has addressed the use of the wavelet decomposition in digital communications. For instance, Wornell (8) has developed the concept of fractal modulation for multiscale communication. Moreover, in works such as (9) and (10), an optimal wavelet decomposition is chosen to account for specific types of channel conditions or transmitter imperfections. A particular problem with much previous work in the area of multiscale communications is that the issue of equalization of multipath channels at the receiver is oftentimes not specifically addressed. This is most likely due to the difficulty in trying to adaptively equalize the channel in two dimensions. Thus it would be desirable to be able to match a particular wavelet to instantaneous channel conditions with a minimal amount of interaction (i.e. feedback) between the receiver and transmitter. However, previous work has not addressed taking advantage of compact realizations of large wavelet families so as to match a wavelet with wireless channel conditions based on the selection of one or more scalar values. In this work, based on well-known compact wavelet decompositions, a parameterized wavelet modulation method is developed in which parameters are chosen to best match the wireless channel conditions.

Wavelet decompositions are normally defined in the continuous domain, where the so-called scaling function φ(x) is first derived (2),(3), φ(x)=Σc _(k)φ(2x−k)   (1) where {c_(k)} is a real sequence. The sequence {c_(k)} is of even length and must satisfy the following: Σc_(k)=2 Σc _(k) c _(k+2m) =2δ( m)   (2)

Another important characteristic of wavelets which determines the “smoothness”, or continuity of the sequence defined by {c_(k)} is the number of vanishing moments. If the wavelet has M (M≧1) vanishing moments, then the following holds: Σ(−1)^(k) k ^(m) c _(k)=0, m=0,1(M−1)   (3)

A corresponding wavelet can now be defined as ψ(x)=Σd _(k)ψ(2x-k)   (4) where d _(k)=(−1)^(k) c _(1-k)   (5)

Thus the dilates and translations of the wavelet function form an orthonormal basis: {{square root}{square root over (2^(j))}ψ(2^(j)x-k)}  (6)

Since the wavelet has compact support, the sequence {c_(k)} is of finite length, assume that the sequence length is 2N. Then the discrete wavelet transform may be defined starting with the two equal length sequences {c_(k)} and {d_(k)}. These two sequences can also be thought of as filters; they collectively form a perfect reconstruction filter bank.

A parameterized construction for wavelet and scaling filters for arbitrary values of N with M≦N vanishing moments was proposed in (3 and (4). Let us assume that for a value of N, the filter coefficients are now denoted as {c_(k) ^(N)}. Given an N-length wavelet parameter set {α_(i)} (−π≦α_(i)<π, 0≦i<N), the coefficients {c_(k) ^(N)} are derived by the recursion $\begin{matrix} \begin{matrix} {c_{0}^{0} = \frac{1}{\sqrt{2}}} \\ {c_{1}^{0} = \frac{1}{\sqrt{2}}} \\ {c_{k}^{n} = {\frac{1}{2}\left\lbrack {{\left( {c_{k - 2}^{n - 1} + c_{k}^{n - 1}} \right)\left( {1 + {\cos\quad\alpha_{n - 1}}} \right)} +} \right.}} \\ \left. {\left( {c_{{2{({n + 1})}} - k - 1}^{n - 1} + c_{{2{({n + 1})}} - k - 3}^{n - 1}} \right)\left( {- 1} \right)^{k}\sin\quad\alpha_{n - 1}} \right\rbrack \end{matrix} & (7) \end{matrix}$

The wavelet construct in (7) is restrictive in the sense that the parameter set cannot in general be defined on the [−π,π]^(N) continuum and still yield a wavelet with at least one vanishing moment. Pollen proved however that wavelets can be defined on the continuum [−π,π]^(N) for arbitrary N (at the cost of smoothness). To examine these types of constructs for a given N, let us define the filter bank matrix F_(N) as $\begin{matrix} {F_{N} = \begin{bmatrix} c_{0} & \ldots & {\quad c_{{2N} - 1}} \\ d_{0} & \ldots & d_{{2N} - 1} \end{bmatrix}} & (8) \end{matrix}$

Then the filter bank matrix for N=1 is $\begin{matrix} {{F_{1} = {\frac{1}{\sqrt{2}}\begin{bmatrix} {\cos\quad\alpha_{0}} & {{- \sin}\quad\alpha_{0}} \\ {\sin\quad\alpha_{0}} & {\cos\quad\alpha_{0}} \end{bmatrix}}}\quad} & (9) \end{matrix}$

The matrix in (9) is also sometimes known as the Givens rotation matrix. Similarly, the filter bank matrix for N=2 is 6) $\begin{matrix} {F_{2} = {\frac{1}{\sqrt{8}}{\quad\begin{bmatrix} {1 - {\cos\quad\alpha_{0}} + {\sin\quad\alpha_{0}}} & {1 + {\cos\quad\alpha_{0}} + {\sin\quad\alpha_{0}}} & {1 + {\cos\quad\alpha_{0}} - {\sin\quad\alpha_{0}}} & {1 - {\cos\quad\alpha_{0}} - {\sin\quad\alpha_{0}}} \\ {{1 - {\cos\quad\alpha_{0}} - {\sin\quad\alpha_{0}}}\quad} & {- \left( {1 + {\cos\quad\alpha_{0}} - {\sin\quad\alpha_{0}}} \right)} & {1 + {\cos\quad\alpha_{0}} + {\sin\quad\alpha_{0}}} & {- \left( {1 - {\cos\quad\alpha_{0}} + {\sin\quad\alpha_{0}}} \right)} \end{bmatrix}}}} & (10) \end{matrix}$

The filter bank expressions for N=2 and N=3 are also sometimes known as Pollen filters, due to the fact that Pollen first proposed these two representations 7). Similarly, the filter bank matrix F₁ becomes the Haar matrix when α₀=π/4 and F₂ is the Daubechies 4-tap filter bank when α₀=π/6. Although the filter selectivity improves with increasing N, this comes at the cost of having to determine a larger set of parameters to define the wavelet. This can be seen in the parameterized expression for {c_(k) ³}: $\begin{matrix} {{c_{0}^{3} = {\frac{1}{4}\left\lbrack {{\left( {1 + {\cos\quad\alpha_{0}} + {\sin\quad\alpha_{0}}} \right)\left( {1 - {\cos\quad\alpha_{1}} - {\sin\quad\alpha_{1}}} \right)} + {2\cos\quad\alpha_{0}\sin\quad\alpha_{1}}} \right\rbrack}}{c_{1}^{3} = {\frac{1}{4}\left\lbrack {{\left( {1 - {\cos\quad\alpha_{0}} + {\sin\quad\alpha_{0}}} \right)\left( {1 + {\cos\quad\alpha_{1}} - {\sin\quad\alpha_{1}}} \right)} - {2\cos\quad\alpha_{0}\sin\quad\alpha_{1}}} \right\rbrack}}{c_{2}^{3} = {\frac{1}{2}\left\lbrack {1 + {\cos\left( {\alpha_{0} - \alpha_{1}} \right)} + {\sin\left( {\alpha_{0} - \alpha_{1}} \right)}} \right\rbrack}}{c_{3}^{3} = {\frac{1}{2}\left\lbrack {1 + {\cos\left( {\alpha_{0} - \alpha_{1}} \right)} - {\sin\left( {\alpha_{0} - \alpha_{1}} \right)}} \right\rbrack}}{c_{4}^{3} = {1 - c_{2}^{3} - c_{0}^{3}}}{c_{5}^{3} = {1 - c_{3}^{3} - c_{1}^{3}}}} & (11) \end{matrix}$

Now two parameters must be determined before specifying the filter bank. In fact, although filter selectivity improves with an increasing number of coefficients, the complexity involved in setting the necessary parameters to form the filter bank also increases.

The wavelet decomposition can now be specified in terms of series of filter banks and resampling stages. Given an input sequence a_(i)(n), then the output sequence may be derived as per the processing depicted in FIG. 2.

In modern technology, the filter process is performed digitally with a computational system such as a general purpose computer or an integrated circuit adapted for digital signal processing.

The number of resampling stages in a wavelet decomposition is sometimes referred to as the number of dilations. This processing can also be represented as a transformation of the input sequence by a unitary matrix. Assume that the input sequence at time index i, a_(i)(n) is of (even) length K (i.e. 0≦n<K) and we want to define a discrete wavelet transformation (DWT) matrix T_(K) of size K by K for a particular filter bank matrix F_(N). Moreover, assume that there are P dilations desired in the transformation. Then using the cited construction, an iterative method for deriving the transformation matrix may be found. Defining the time scale index as 1(0≦1≦P), a K-by-K filter bank matrix can be defined for each time scale: $\begin{matrix} {{C_{l} = \begin{bmatrix} {\left. \begin{matrix} F_{N} & \quad & 0_{2{x{({K - {2N} + 1})}}} \\ 0_{2{x2}} & {\overset{.}{F}}_{N} & 0_{2{x{({K - {2N} - 1})}}} \\ 0_{2{x4}} & F_{N} & 0_{2{x{({K - {2N} - 3})}}} \\ \vdots & ⋰ & \vdots \end{matrix} \right\}\left( {K - R} \right){rows}} \\ \begin{matrix} 0_{RxR} & \quad & \quad & \quad & \quad & I_{R} & \quad & \quad \end{matrix} \end{bmatrix}}{R = \left\{ \begin{matrix} 0 & {l = 0} \\ {\sum\limits_{i = 1}^{P}\frac{K}{2^{i}}} & {l > 0} \end{matrix} \right.}} & (12) \end{matrix}$

In (12), 0_(m×n) is the zero matrix of m rows by n columns, and I_(R) is the identity matrix of R rows by R columns. For each dilation, a permutation matrix P_(υ)(1≦υ≦P) can be defined as well: $\begin{matrix} {{P_{v} = \begin{bmatrix} {\left. \begin{matrix} 1 & 0 & \ldots & 0 & 0 & 0 & 0 & 0 \\ 0 & 0 & 1 & 0 & \ldots & 0 & 0 & 0 \\ 0 & 0 & 0 & 0 & 1 & 0 & \ldots & 0 \\ \vdots & \quad & \vdots & \quad & \quad & \vdots & \quad & \vdots \end{matrix} \right\}\left( \frac{K - R}{2} \right){rows}} \\ {\left. \begin{matrix} 0 & 1 & 0 & \ldots & 0 & 0 & 0 & 0 \\ 0 & 0 & 0 & 1 & 0 & \ldots & 0 & 0 \\ 0 & 0 & 0 & 0 & 0 & 1 & \ldots & 0 \\ \vdots & \quad & \vdots & \quad & \quad & \vdots & \quad & \vdots \end{matrix} \right\}\left( \frac{K - R}{2} \right){rows}} \\ \begin{matrix} 0_{RxR} & \quad & I_{R} & \quad & \quad & \quad \end{matrix} \end{bmatrix}}{R = \left\{ \begin{matrix} 0 & {l = 0} \\ {\sum\limits_{i = 1}^{P}\frac{K}{2^{i}}} & {l > 0} \end{matrix} \right.}} & (13) \end{matrix}$

Thus, for P dilations, the unitary transform matrix T_(K) (P) may be determined as T _(K)(P)=C _(P) P _(P) . . . C ₁ C ₁ C ₀   (14)

This matrix may now be used to modulate an input vector a_(i)=[a_(i)(0) . . . a_(i)(K-1)]^(T) (i being the symbol index). Illustratively, the elements of the input signal vector ai are the coefficients on a set of basis functions that represent the data, e.g. speech.

Thus the output sequence of such a modulation can be formed as the result of a matrix-vector multiplication x_(i)=T_(K) ^(T)(P)a_(i). As described above, the single net matrix T_(K)(P) represents P inner products of the wavelet with the data. If we assume that F_(N) results from a compact realization of a wavelet, then we can also assume that T_(K)(P) is a function of the wavelet parameter set {α_(i)}, i.e. x_(i)=T_(K)(P,{α_(i)})a_(i).

SUMMARY OF THE INVENTION

The invention relates to a multi-carrier communications system that adaptively selects a set of wavelets that match channel conditions.

A feature of the invention is the use of a compact parameterization that permits the generation of an infinite number of wavelet and scaling filter pairs using a finite set of parameters such as the number of dilations and the filter parameters.

Another feature of the invention is the elimination of an equalizer in the receiver by altering the parameters of the wavelets to compensate in advance for changing channel conditions.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 shows time-frequency relationship in a wave packet system.

FIG. 2 shows a block diagram of a wavelet transmission system.

FIG. 3 shows quantization noise characteristics

FIG. 4 shows feedback effects in a 2-tap wavelet system.

FIG. 5 shows feedback effects in a 4-tap wavelet system.

FIG. 6 shows dilation effects in a 2-tap wavelet system.

FIG. 7 shows dilation effects in a 4-tap wavelet system.

FIG. 8 compares BER between a wavelet and an OFDM system.

FIG. 9 compares BER in fixed versus adaptive wavelet systems, 2-tap.

FIG. 10 compares BER in fixed verses adaptive wavelet systems, 4-tap.

FIG. 11 shows BER in various systems.

DETAILED DESCRIPTION

A wavelet-based transmission system can be formed starting with an input vector a_(i)=[a_(i)(0) . . . a_(i)(K-1)]^(T) and forming an output vector x_(i)=T_(K) ^(T)(P,{α_(i)})a_(i) in a similar fashion to an OFDM system. However, unlike OFDM systems, which can employ simple equalization structures through the use of cyclic convolution properties of the underlying discrete Fourier transform modulation, a wavelet-based system can be sensitive to frequency-selective wireless channels and therefore the equalization problem for wavelets can become complex.

The output vector x_(i) is transmitted serially. If it is assumed that the wireless channel can be described by an L-tap channel vector h_(i)=[h_(i)(0) . . . h_(i)(L-1)]^(T) (h_(i) is assumed to have unit norm), then the elements of the received signal y_(i) may be represented as (assuming L<K) $\begin{matrix} \begin{matrix} {{y_{i}(k)} = {{\sum\limits_{m = 0}^{k}{{h_{i}(m)}x_{i}\left( {k - m} \right)}} + {\sum\limits_{m = {k + 1}}^{L - 1}{h_{i}(m){x_{i - 1}\left( {K + k - m} \right)}}} +}} \\ {{n_{i}(k)},{0 \leq k \leq {L - 1}}} \\ {{= {{\sum\limits_{m = 0}^{L - 1}{{h_{i}(m)}{x_{i}\left( {k - m} \right)}}} + {n_{i}(k)}}},{L \leq k < {K - 1}}} \end{matrix} & (15) \end{matrix}$

In (15), n_(i)(k) is an additive noise term. The typical equalization techniques for this type of received signal are decision-aided, requiring an estimate of the previous transmitted symbol a_(i-1) to account for intersymbol interference (ISI) and estimates of the individual entries of x_(i) to account for intercarrier interference (ICI). This expression can also be expressed in matrix-vector format: y _(i) =H _(i) x _(i) +H _(ISI) _(—) _(i) x _(i-1) +n _(i)   (16)

In (16), H_(i), H_(ISI) _(—) _(i), and n_(i) are described as below: $H_{i} = \begin{bmatrix} h_{0} & 0 & \ldots & 0 & 0 & 0 & 0 & 0 & 0 \\ h_{1} & h_{0} & 0 & {\ldots~} & 0 & 0 & 0 & 0 & 0 \\ h_{2} & h_{1} & h_{0} & 0 & \ldots & 0 & 0 & 0 & 0 \\ \vdots & \quad & \quad & \quad & \vdots & \quad & \quad & \quad & \vdots \\ 0 & \ldots & 0 & 0 & 0 & h_{L - 1} & \ldots & h_{1} & h_{0} \end{bmatrix}$ $H_{ISI\_ i} = \begin{bmatrix} 0 & \ldots & 0 & h_{L - 1} & \ldots & h_{2} & h_{1} \\ 0 & 0 & \ldots & 0 & h_{L - 1} & \ldots & h_{2} \\ 0 & 0 & 0 & \ldots & 0 & ⋰ & \vdots \\ 0 & 0 & 0 & 0 & \ldots & 0 & h_{L - 1} \\ 0 & 0 & 0 & 0 & 0 & \ldots & 0 \\ \vdots & \quad & \quad & \vdots & \quad & \quad & \vdots \\ 0 & 0 & 0 & 0 & 0 & \ldots & 0 \end{bmatrix}$ $n_{i} = \begin{bmatrix} {n_{i}(0)} & {n_{i}(1)} & \ldots & {n_{i}\left( {K - 1} \right)} \end{bmatrix}^{T}$

What is proposed instead is to utilize compact parameterization of wavelets to pre-equalize the channel. In other words, if there is an estimate of a set of parameters for a given filter bank matrix F_(N) which maximizes received signal quality, then this information can be used to modify the modulation matrix T_(K). If it is assumed that there is sufficient training data to form a channel estimate at the receiver ĥ_(i), then the best wavelet parameterization set may be found utilizing additional training data and the K-by-K channel estimation matrix: $\begin{matrix} {{\hat{H}}_{i} = \begin{bmatrix} {\hat{h}}_{0} & 0 & \ldots & 0 & 0 & 0 & 0 & 0 & 0 \\ {\hat{h}}_{1} & {\hat{h}}_{0} & 0 & {\ldots~} & 0 & 0 & 0 & 0 & 0 \\ {\hat{h}}_{2} & {\hat{h}}_{1} & {\hat{h}}_{0} & 0 & \ldots & 0 & 0 & 0 & 0 \\ \vdots & \quad & \quad & \quad & \vdots & \quad & \quad & \quad & \vdots \\ 0 & \ldots & 0 & 0 & 0 & {\hat{h}}_{L - 1} & \ldots & {\hat{h}}_{1} & {\hat{h}}_{0} \end{bmatrix}} & (17) \end{matrix}$

Clearly, applying a unitary transformation will not reduce the power of additive Gaussian noise; however, a unitary transformation such as a DWT could be applied to minimize the two sources of interference seen in frequency selective channels by multicarrier systems, namely ICI and ISI. Using this channel estimation matrix, the optimal wavelet parameterization set which minimizes ICI for a given P can be found as $\begin{matrix} {\min\limits_{{\{\alpha_{r}\}} \in \mathcal{R}^{N}}{{{{T_{K}^{T}\left( {P,\left\{ \alpha_{r} \right\}} \right)}{\hat{H}}_{i}{T_{K}\left( {P,\left\{ \alpha_{r} \right\}} \right)}} - I_{K}}}} & (18) \end{matrix}$

Additional criteria could be established to minimize ISI. This would involve selection of a wavelet transform that would render ISI contributions from any given wavelet symbol to the ensuing wavelet symbol negligible. Based on the ISI contributions evident in (15), we define the K-by-K channel matrix: $\begin{matrix} {{\hat{H}}_{ISI\_ i} = \begin{bmatrix} 0 & \ldots & 0 & {\hat{h}}_{L - 1} & \ldots & {\hat{h}}_{2} & {\hat{h}}_{1} \\ 0 & 0 & \ldots & 0 & {\hat{h}}_{L - 1} & \ldots & {\hat{h}}_{2} \\ 0 & 0 & 0 & \ldots & 0 & ⋰ & \vdots \\ 0 & 0 & 0 & 0 & \ldots & 0 & {\hat{h}}_{L - 1} \\ 0 & 0 & 0 & 0 & 0 & \ldots & 0 \\ \vdots & \quad & \quad & \vdots & \quad & \quad & \vdots \\ 0 & 0 & 0 & 0 & 0 & \ldots & 0 \end{bmatrix}} & (19) \end{matrix}$

Thus, another optimization criterion that minimizes ISI may be considered: $\begin{matrix} {\min\limits_{{\{\alpha_{r}\}} \in \mathcal{R}^{N}}{{{T_{K}^{T}\left( {P,\left\{ \alpha_{r} \right\}} \right)}{\hat{H}}_{ISI\_ i}{T_{K}\left( {P,\left\{ \alpha_{r} \right\}} \right)}}}} & (20) \end{matrix}$

In addition to selection of a wavelet that minimizes ICI and ISI, the wavelet itself cannot be considered to have approximately equalized the channel unless the energy along the resultant diagonal is maximized. This results directly in maximizing received signal energy. The criterion that maximizes the energy along the diagonal (assuming Λ is a diagonal operator) is: $\begin{matrix} {\min\limits_{{\{\alpha_{r}\}} \in \mathcal{R}^{N}}\frac{1}{{\Lambda\left( {{T_{K}^{T}\left( {P,\left\{ \alpha_{r} \right\}} \right)}{\hat{H}}_{i}{T_{K}\left( {P,\left\{ \alpha_{r} \right\}} \right)}} \right)}}} & (21) \end{matrix}$

It cannot be assumed that a single parameter set will minimize ISI and ICI while maximizing recovered channel energy. Thus, the wavelet parameter selection should be based on a criterion that minimizes all of the residual interference: $\begin{matrix} {{\min\limits_{{\{\alpha_{r}\}} \in \mathcal{R}^{N}}{{{{T_{K}^{T}\left( {P,\left\{ \alpha_{r} \right\}} \right)}{\hat{H}}_{i}{T_{K}\left( {P,\left\{ \alpha_{r} \right\}} \right)}} - {{\hat{h}}_{0}I_{K}{ + }{T_{K}^{T}\left( {P,\left\{ \alpha_{r} \right\}} \right)}{\hat{H}}_{ISI\_ i}{T_{K}\left( {P,\left\{ \alpha_{r} \right\}} \right)}}}}} + \frac{1}{{\Lambda\left( {{T_{K}^{T}\left( {P,\left\{ \alpha_{r} \right\}} \right)}{\hat{H}}_{i}{T_{K}\left( {P,\left\{ \alpha_{r} \right\}} \right)}} \right)}}} & (22) \end{matrix}$

Finally, the wavelet transmission method is as depicted in FIG. 2.

Summarizing, the process according to the invention is:

Starting with an estimated channel matrix, and an initial parameter set alpha that specifies an initial wavelet (assuming P is fixed).

Send training signals from the transmitter (base station) to the receiver.

In the receiver iterate (or otherwise compute) the value of adjusted wavelet parameters that minimize the ICI, the ISI and the total residual interference.

Send the adjusted parameters back along a feedback path to the transmitter.

Use the adjusted parameters in transmission for the next period until the next adjustment.

It should be noted that there could still be residual interference even after wavelet selection. This would result from effects such as feedback latency with respect to the coherence time of the channel, reduced parameter search space, etc. As a result, some limited form of interference cancellation may still be necessary even when the wavelet is best matched to channel conditions. Moreover, using the metric in (22) does not maximize diversity seen in the system; it only optimizes ICI and ISI.

The design of parameterized wavelets poses a classic problem in adaptive waveform design for wireless transmission in a system where the receiver must convey information on the channel conditions to the transmitter: how to minimize the throughput required for relaying the necessary information from the receiver to the transmitter for accurate waveform selection. This is normally not a problem in time division duplexed (TDD) systems, as it is generally assumed that a TDD transmitter can estimate channel conditions seen at the TDD receiver without feedback. However, paired-band systems generally do not have sufficient correlation between transmit and receive frequencies used by a given transceiver. Thus minimizing the amount of information to be relayed by the receiver to the transmitter for waveform selection is critical.

However, the wavelet filter bank matrix selectivity clearly improves with increasing order 2N. This requires greater parameterization, thus potentially increasing the required feedback information. Thus given a maximum feedback data payload R (bits), there exists a tradeoff between increasing the filter bank selectivity and reducing the quantization noise of the feedback parameters. If we assume that each value in the parameter set {α_(i)} is uniformly quantized between 0 and π, we can find the mean quantization error as $\begin{matrix} {\sigma_{q}^{2} = {\frac{1}{12}\left( \frac{\pi}{2^{\lfloor\frac{R}{N}\rfloor}} \right)^{2}}} & (23) \end{matrix}$

In (23), the quantization error is derived from the classic result on uniform quantization of a uniformly distributed random variable. Unfortunately, this error variance in and of itself may not give sufficient insight into the degradation of a wavelet-based communications system due to throughput limitations for feedback of wavelet parameter selection. However, it could provide an indication of the optimal filter bank order given specific feedback limitations. The quantization noise for filter orders of 2 and 4 is depicted in FIG. 4. In the left region of the curve, quantization noise increases to unacceptable levels with increasing filter order. However, in the right region of the curve, quantization noise tends to drop below unacceptable levels with increasing filter order (although quantization noise still worsens with increasing filter order for any given value of R). This is not surprising, as it stands to reason that increasing throughput payloads will allow for more accurate quantization of an increasing number of wavelet parameters, ergo allowing for higher-order filterbanks.

A relevant question is what constitutes an “unacceptable level” of quantization noise, as it would stand to reason that for a given value of R, the filter order that should be selected is the maximum one which falls below this level.

FIG. 2 illustrates the process of forming the coefficients that produce a given wavelet. A set of signal coefficients a_(i) representing the speech or other data are input on the left and are processed (e.g. by a general purpose computer system) to generate the final coefficients that define the particular wavelets used in the transmitter.

In operation, in a system according to the invention, the transmitter periodically sends a reference signal (a training sequence) to the receiver. The receiver applies known algorithms to the received signal to estimate the channel matrix. The receiver then transforms the channel matrix using various trial parameters and selects the parameters that give the result that best fits a criterion such as that expressed in Equation 22 that minimizes ICI and ISI.

The “best fit” parameters are relayed to the transmitter along a feedback channel having a limited number of bits

The transmitter then prepares packets using the parameters sent to it by the receiver for the duration of the period.

A test over a specific channel profile was conducted for several different types of wavelet decompositions to examine the feasibility of adapting the wavelet to the instantaneous channel conditions based on their parameterization. Assuming an input vector of dimension 32 modulated by a square wavelet transformation matrix, a transmission system was tested at a wavelet symbol rate of 250 kHz (symbol rate in this case refers to the rate at which all 32 wavelet coefficients resulting from a single input vector was transmitted). The two-tap and four-tap wavelet models of (9) and (10) were examined for performance of BPSK signaling in a multipath fading channel (channel tap power profile of [0.8 0.1 0.1]) at an assumed carrier frequency of 5 GHz and velocity of 3 km/hr. In the system modeling, the wavelet parameter was selected every 50 symbols. Under these conditions, the performance of the 2-tap and 4-tap wavelets with respect to feedback quantization is given in FIG. 5 and FIG. 6, respectively.

Under these conditions, the benefits of increasing feedback naturally become less prevalent as the number of bits increases. In this environment, although increasing feedback would yield 6 dB of quantization noise relief, this clearly does not translate to the same benefits in terms of overall performance. It should also be noted that the results are graphed in terms of raw BPSK bit error rate; an additional error correcting code could translate to less benefits for increasing the number of feedback bits. Another thing to notice is that the performance of the 2-tap and 4-tap wavelet systems was roughly the same.

In addition, the number of dilations could possibly impact the performance of the proposed method of wavelet transmission. Once again, the results for both the 2-tap and 4-tap wavelets assuming 3-bit feedback quantization are provided in FIG. 7 and FIG. 8, respectively. Increasing the number of dilations did not provide much benefit in improving performance of the system. This could most likely be because of the fact that the channel conditions under consideration are already fairly compact in frequency and in time; the greatest performance enhancement of this method is to change the fundamental wavelet filterbank rather than changing the time-frequency decomposition of the transmitted signal.

To better appreciate the maximum achievable benefits of multiscale communications, an OFDM system was simulated using identical channel conditions and identical input vector size. A comparison to the 2-tap wavelet with and one bit of feedback is given in FIG. 9. The performance enhancements due to wavelet modulation are seen to be as much as 3 dB. Taking into account that the OFDM system is transmitted in quadrature while the wavelet system does not require quadrature transmission for BPSK signaling, it is clear that the spectral efficiency of the wavelet system is potentially much higher than the OFDM system.

Some more results were obtained for the adaptive method in fading channels and compared to the use of a fixed method (in this case, the Haar basis and Daubechies 4-tap). Using a channel tap profile of [0.8 0.2] along with a symbol rate of 125 kHz and an input vector size of 64, the fixed and adaptive methods were compared under fading conditions with a carrier frequency of 5 GHz and a velocity of 150 km/hr. A rate ½ convolutional code was used, meaning that an input segment of 32 bits was coded into an input vector of 64 binary symbols (BPSK signaling). The results are shown in FIG. 10 and FIG. 11.

Both of the methods seem to converge at low SNR. This is due to AWGN being the overriding source of error in this operating region; no adaptation can enhance the link under such conditions. In addition, note that both methods do result in an error floor. For the fixed case, this is due to a suboptimal wavelet being used for transmission. In the adaptive case, this is due in part to the use of limited feedback. However, even assuming infinitely precise quantization of the parameterization space, it should still be re-asserted that the selected wavelet will be an approximation of the equalizing transform.

The wavelet selection method given above, particularly the criterion presented in Equation (22), was examined versus the fixed method in the static channel profile of [0.5 0.3 0.1 0.1]. In this case, the feedback quantization was increased to 5 bits so as to more accurately examine whether the wavelet selection actually selects the best wavelet. More specifically, under the assumptions of coarse quantization, the likelihood of picking the correct wavelet using the wavelet selection method increases; therefore the accuracy of this method is not so easy to assess under such conditions. The results for the 2-tap method are given in FIG. 12. Note that the wavelet selection method retained the gains of the adaptive method versus the use of a fixed wavelet. Note also that the error floor is still present, although at a bit error rate between 10⁻⁴ and 10⁻³.

Although the invention has been described with respect to a limited number of embodiments, those skilled in the art will appreciate that other embodiments may be constructed within the spirit and scope of the following claims.

-   (1) Vetterli, Martin and Jelena Kovacevic. Wavelets and Subband     Coding. Englewood Cliffs, N J: Prentice-Hall Inc., 1995. -   (2) Strang, Gilbert and Truong Nguyen. Wavelets and Filter Banks.     Wellesley, Mass.: Wellesley-Cambridge Press, 1996. -   (3) Zou, Hehong and Ahmed H. Tewfik. “Parameterization of Compactly     Supported Orthonormal Wavelets.” IEEE Transactions on Signal     Processing. Vol. 41. No. 3. March 1993. pp. 1428-1431. -   (4) Schneid, J. and S. Pittner. “On the Parameterization of the     Coefficients of Dilation Equations for Compactly Supported     Wavelets.” Computing. Vol.51. May 1993. pp.165-173. -   (5) Vidakovic, Brani. “Pollen Bases and Daubechies-Lagarias     Algorithm in MATLAB.” Jacket's Wavelets website.     http://www.isye.gatech.edu/˜brani/datasoft/DL.pdf. -   (6) Silva, Vitor and Luis de Sa. “Analytical Optimization of CQF     Filter Banks.” IEEE Transactions on Signal Processing. Vol. 44.     No. 6. June 1996. pp. 1564-1568. -   (7) Pollen, D. “Parameterization of Compactly Supported Wavelets.”     Aware Inc. technical report AD890503. May 1989. -   (8) Wornell, Gregory W. “Emerging Applications of Multirate Signal     Processing and Wavelets in Digital Communications.” Proceedings of     the IEEE. Vol. 84. No. 4. April 1996. pp. 586-603. -   (9) Lindsey, Alan R. “Wavelet Packet Modulation for Orthogonally     Multiplexed Communication.” IEEE Transactions on Signal Processing.     Vol. 45. No. 5. May 1997. pp. 1336-1339. -   (10) Wong, K.-M., Wu, J., Davidson, T. N. and Jin, Q “Wavelet Packet     Division Multiplexing and Wavelet Packet Design under Timing Error     Effects.” IEEE Transactions on Signal Processing. Vol. 45. No. 12.     December 1997. pp. 2877-2890. 

1. A telecommunications system comprising: means for receiving an input vector signal; controllable variable transformation means for modulating said input signal with a transform having a compact parameter representation with P dilations, and further characterized by a parameter set, said transformations means having a parallel output set of data; parallel to serial means for converting said parallel set of data to a serial set of data; and transmitting means for transmitting said serial set of data along a channel; receiving means for receiving said serial set of data; serial to parallel means for converting received serial data to parallel data; and inverse transformation means for demodulating said parallel data in accordance with said compact parameter representation.
 2. A system according to claim 1, in which said process of modulating is effected by multiplying a vector representing input data with a modulation matrix representing the transform.
 3. A system according to claim 1, in which said receiving means contains means for processing test data to extract parameters for specifying said modulation matrix.
 4. A system according to claim 3, in which said process of modulating is effected by multiplying a vector representing input data with a modulation matrix representing the transform.
 5. A system according to claim 3, in which said means for processing test data to extract parameters seeks parameters that minimize residual interference, including ISI and ISI.
 6. A system according to claim 3, in which said means for processing test data to extract parameters seeks parameters that minimize ICI.
 7. A system according to claim 3, in which said means for processing test data to extract parameters seeks parameters that minimize ISI.
 8. A telecommunications system comprising: means for receiving an input vector signal; controllable variable transformation means for modulating said input signal with a transform having a compact parameter representation with P dilations, and further characterized by a parameter set, said transformations means having a parallel output set of data; parallel to serial means for converting said parallel set of data to a serial set of data; transmitting means for transmitting said serial set of data along a channel; receiving means for receiving said serial set of data; serial to parallel means for converting received serial data to parallel data; inverse transformation means for demodulating said parallel data in accordance with said compact parameter representation; and adjustment means for adjusting said parameter set in response to changing channel conditions.
 9. A system according to claim 8, in which said process of modulating is effected by multiplying a vector representing input data with a modulation matrix representing the transform.
 10. A system according to claim 8, in which said receiving means contains means for processing test data to extract parameters for specifying said modulation matrix.
 11. A system according to claim 10, in which said process of modulating is effected by multiplying a vector representing input data with a modulation matrix representing the transform.
 12. A system according to claim 10, in which said means for processing test data to extract parameters seeks parameters that minimize residual interference, including ISI and ISI.
 13. A system according to claim 10, in which said means for processing test data to extract parameters seeks parameters that minimize ICI.
 14. A system according to claim 10, in which said means for processing test data to extract parameters seeks parameters that minimize ISI.
 15. A method of operating a telecommunications system comprising the steps of: receiving an input vector signal; modulating said input signal in controllable variable transformation means with a transform having a compact parameter representation with P dilations, and further characterized by a parameter set, said transformation means having a parallel output set of data; converting said parallel set of data to a serial set of data; transmitting said serial set of data along a channel; receiving said serial set of data; converting received serial data to parallel data; and demodulating said parallel data in accordance with said compact parameter representation.
 16. A method according to claim 15, in which said receiving means processes test data to extract parameters for specifying said modulation matrix.
 17. A method according to claim 15, in which said process of modulating is effected by multiplying a vector representing input data with a modulation matrix representing the transform.
 18. A method according to claim 15, in which said step of processing test data to extract parameters seeks parameters that minimize residual interference, including ISI and ISI.
 19. A method according to claim 15, in which said step of processing test data to extract parameters seeks parameters that minimize ICI.
 20. A method according to claim 15, in which said step of processing test data to extract parameters seeks parameters that minimize ISI.
 21. A method of operating a telecommunications system comprising the steps of: receiving an input vector signal; modulating said input signal in controllable variable transformation means with a transform having a compact parameter representation with P dilations, and further characterized by a parameter set, said transformation means having a parallel output set of data; converting said parallel set of data to a serial set of data; adjusting said parameter set in response to changing channel conditions; transmitting said serial set of data along a channel; receiving said serial set of data; converting received serial data to parallel data; and demodulating said parallel data in accordance with said compact parameter representation.
 22. A method according to claim 21, in which said receiving means processes test data to extract parameters for specifying said modulation matrix.
 23. A method according to claim 21, in which said process of modulating is effected by multiplying a vector representing input data with a modulation matrix representing the transform.
 24. A method according to claim 21, in which said step of processing test data to extract parameters seeks parameters that minimize residual interference, including ISI and ISI.
 25. A method according to claim 21, in which said step of processing test data to extract parameters seeks parameters that minimize ICI.
 26. A method according to claim 21, in which said step of processing test data to extract parameters seeks parameters that minimize ISI.
 27. An article of manufacture comprising a program storage medium readable by a computer, the medium embodying instructions executable by the computer for performing method steps of operating a telecommunications system comprising the steps of: receiving an input vector signal; modulating said input signal in controllable variable transformation means with a transform having a compact parameter representation with P dilations, and further characterized by a parameter set, said transformation means having a parallel output set of data; converting said parallel set of data to a serial set of data; transmitting said serial set of data along a channel; receiving said serial set of data; converting received serial data to parallel data; and demodulating said parallel data in accordance with said compact parameter representation.
 28. A method according to claim 15, in which said receiving means processes test data to extract parameters for specifying said modulation matrix.
 29. A method according to claim 15, in which said process of modulating is effected by multiplying a vector representing input data with a modulation matrix representing the transform.
 30. A method according to claim 15, in which said step of processing test data to extract parameters seeks parameters that minimize residual interference, including ISI and ISI.
 31. A method according to claim 15, in which said step of processing test data to extract parameters seeks parameters that minimize ICI.
 32. A method according to claim 15, in which said step of processing test data to extract parameters seeks parameters that minimize ISI.
 33. An article of manufacture comprising a program storage medium readable by a computer, the medium embodying instructions executable by the computer for performing method steps of operating a telecommunications system comprising the steps of: receiving an input vector signal; modulating said input signal in controllable variable transformation means with a transform having a compact parameter representation with P dilations, and further characterized by a parameter set, said transformation means having a parallel output set of data; converting said parallel set of data to a serial set of data; adjusting said parameter set in response to changing channel conditions; transmitting said serial set of data along a channel; receiving said serial set of data; converting received serial data to parallel data; and demodulating said parallel data in accordance with said compact parameter representation.
 34. A method according to claim 33, in which said receiving means processes test data to extract parameters for specifying said modulation matrix.
 35. A method according to claim 33, in which said process of modulating is effected by multiplying a vector representing input data with a modulation matrix representing the transform.
 36. A method according to claim 33, in which said step of processing test data to extract parameters seeks parameters that minimize residual interference, including ISI and ISI.
 37. A method according to claim 33, in which said step of processing test data to extract parameters seeks parameters that minimize ICI.
 38. A method according to claim 33, in which said step of processing test data to extract parameters seeks parameters that minimize ISI.
 39. A transmitter comprising: an input port for receiving an input vector signal to be transmitted; a modulator for modulating said input signal with a transform having a compact parameter representation with P dilations, and further characterized by a parameter set, said modulator providing a parallel output set of data; a parallel to serial converter for converting said parallel set of data to a serial set of data; and a transmit antenna for transmitting said serial set of data along a channel.
 40. A receiver comprising: a receive antenna for receiving a serial set of data over a channel; a serial to parallel converter for converting said received serial data to parallel data; and a demodulator for demodulating said parallel data in accordance with a compact parameter representation having P dilations, said representation further characterized by a parameter set.
 41. The receiver of claim 40 further comprising: a controller for adjusting said parameter set in response to changing conditions of the channel. 